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Abstract. In this work a new theory of nucleation is developed and used to investigate the effect 
of water activity on the formation of ice within supercooled droplets. The new theory is based on 
a novel concept where the interface is assumed to be made of liquid molecules “trapped” by the 
solid matrix. Using this concept new expressions are developed for the critical ice germ size and 
5 the nucleation work, with explicit dependencies on temperature and water activity. However unlike 
previous approaches, the new theory does not depend on the interfacial tension between liquid and 
ice. Comparison against experimental results shows that the new theory is able to reproduce the 
observed effect of water activity on nucleation rate and freezing temperature. It allows for the first 
time a theoretical derivation of the constant shift in water activity between melting and nucleation. 
10 The new theory offers a consistent thermodynamic view of ice nucleation, simple enough to be 
applied in atmospheric models of cloud formation. 

1 Introduction 

Ice formation by freezing of supercoooled droplets is an important natural and technological process. 
In the atmosphere it leads to the formation of cirrus and determines the freezing level of convective 
1 5 clouds (Pruppacher and Klett, 1997). At temperatures below 238 K and in the absence of ice forming 
nuclei, freezing proceeds by homogeneous nucleation. A significant fraction of cirrus in the upper 
troposphere forms by this mechanism (Gettelman et al., 2012; Barahona et al., 2013). Cirrus clouds 
impact the radiative balance of the upper troposphere (Fu, 1996) and play a role in the transport of 
water vapor to the lower stratosphere (e.g., Barahona and Nenes, 2011; Jensen and Pfister, 2004; 
20 Hartmann et al., 2001). Correct parameterization of ice formation is therefore crucial for reliable 
climate and weather prediction (Lohmann and Feichter, 2005). Many experimental and theoretical 
studies have been devoted to the study of homogeneous nucleation (e.g., Kashchiev, 2000; Murray 
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et al., 2010b; Wu et al., 2004, and references therein), yet there are still significant gaps in the 
understanding of ice formation within supercooled droplets. 

The theoretical study of homogeneous ice nucleation is commonly approached by using the Clas- 
sical Nucleation Theory (CNT) (e.g., Khvorostyanov and Curry, 2004) or by molecular dynamics 
(MD) simulations (e.g., Matsumoto et al., 2002). Density functional theory and direct kinetic models 
have also been employed (Laaksonen et al., 1995). MD and other detailed approaches offer a unique 
look to the microscopic mechanism of ice nucleation. However, for climate simulations simplified 
and efficient descriptions of ice nucleation are required. Thus CNT is typically used to generate 
ice formation parameterizations in atmospheric models (Liu and Penner, 2005; Khvorostyanov and 
Curry, 2004). 

CNT is often criticized due to the usage of the so-called capillary approximation, i.e., the assump- 
tion that the properties of ice clusters at nucleation are the same as those of the bulk (Kashchiev, 
2000). Among the parameters used by CNT, the ice-liquid interfacial tension (also called specific 
surface energy), ( 7 ; w , is generally the most uncertain. CNT predictions are highly sensitive to cr iw , 
however direct measurements of a- lw are typically difficult and are only available at low supercooling 
(Pruppacher and Klett, 1997). Furthermore, other factors not considered by CNT like crystal shape, 
type and size, and the lack of a well-defined ice liquid interface may also play a role in determining 
er iw (Wu et al., 2004; MacKenzie, 1997; Kashchiev, 2000; Murray et al., 2010a). Thus er iw is often 
found by fitting to CNT predictions to experimental measurements of the nucleation rate (Murray 
et al., 2010a; Khvorostyanov and Curry, 2004). The caveat is that cr iw obtained by this method often 
differs significantly from theoretical estimates (MacKenzie, 1997), casting doubt into the validity of 
CNT. 

Due to the shortcomings of CNT, experimental data are most often used to describe homogeneous 
freezing in atmospheric models (e.g., Barahona and Nenes, 2008; Karcher and Lohmann, 2002). 
Experimental studies generally agree on the freezing temperature of pure water with typical variation 
below one degree (Murray et al., 2010a; Pruppacher and Klett, 1997). For liquid solutions this 
picture is somewhat complicated by the chemical heterogeneity of liquid droplets in the atmosphere. 
Empirical correlations were often developed based on (NH^SCL and H 2 SO 4 model solutions (e.g., 
Tabazadeh et al., 1997; Jensen et al., 1991). This issue was resolved by Koop et al. (2000) who 
demonstrated that when parameterized in terms of the water activity, a w , freezing temperatures are 
independent of the nature of the solute. Furthermore, the authors showed that when plotted in a 
T — a w diagram, the melting and nucleation curves can be translated by a constant shift in water 
activity. This finding has been confirmed by several independent studies (Marcolli et al., 2007; 
Wang and Knopf, 201 1) and has been be referred to as the “water activity criteria”. The Koop et al. 
(2000) (hereafter K00) parameterization has been incorporated in several global atmospheric models 
(e.g., Barahona et al., 2010; Liu et al., 2007). 

The study of Koop et al. (2000) suggested that a general thermodynamic formulation of ice nucle- 
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ation could be achieved, however such theory has been so far elusive. Current formulations of CNT 
carry a dependency on a w and it has been suggested that CNT can explain the water activity criteria 
(e.g., Khvorostyanov and Curry, 2004). However in these studies er iw is typically modified to adjust 
CNT to K00, therefore the results obtained in this way are not independent of Koop et al. (2000) 
data. In fact, Koop et al. (2000) suggested that CNT and K00 can be reconciled if erj w is allowed to 
vary with a w , although no theoretical support was provided for it. Baker and Baker (2004) took an 
alternative approach and showed that the freezing temperatures measured by K00 were consistent 
with the point of maximum compressibility of water. The authors derived an empirical relation be- 
tween a w and the osmotic pressure which was then to used to determine freezing temperatures. The 
work of Baker and Baker (2004) demonstrated that the water activity criteria can be understood in 
terms of water compressibility as long as certain empirical criteria are met. However a theoretical 
basis for this behavior was not provided. Furthermore, the liquid-only approach used by the authors 
is in contrast with MD simulations showing that ice nucleation originates in differentiated, ice-like 
regions within the liquid phase (Matsumoto et al., 2002). 

In this work a new approach is proposed to describe ice formation by homogeneous nucleation. 
The new model relies on a novel picture of the solid-liquid transition placing emphasis on the entropy 
changes across the interface. The new model is used to analyze the effect of water activity on ice 
formation and ice nucleation rates. 

2 Theory 

Consider the system of Fig. 1 . The liquid droplet is assumed to be large enough so that nucleation is 
more likely to occur within the bulk of the liquid than at the droplet surface. The liquid is assumed to 
be homogeneously mixed and its cluster distribution in steady state. For simplicity it is assumed that 
only two components are present in solution, water (subscript, “ w ”) and a solute (subscript, “y”), 
although this assumption can be easily relaxed if more than one solute is present. The Gibbs free 
energy of the system in stage 1 (before germ formation) is given by 

Nwflw'X ( 1 ) 

where N w and N y are the total number of water and solute molecules present in the liquid phase, 
respectively, and y, w \ and y, yt \ their respective chemical potentials. 

After the formation of the germ (stage 2, Fig. 1) it is advantageous to consider the solid-liquid 
interface as a phase distinct from the bulk (Gibbs et al., 1928). It is assumed that no atoms of y 
are present in the bulk of the solid phase although they may be present at the interface. However 
the dividing surface is selected so that the molecular excess of solute at the interface is zero (this 
is further analyzed in Section 2.1). The assumption of a solute-free solid is justified on molecular 
dynamics simulations showing a rejection of ions into a unfrozen layer of brine away from the germ 
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(Bauerecker et al., 2008). With this, the Gibbs free energy of the system in stage 2 is given by 


G 2 — ( X w Tl s 1^ls)t^w,2 “b AyMy,2 ”b ^sk’w^s T TllsHw,ls (X) 

where n s and ni s are the number of atoms in the bulk of the germ and in the interface, respectively, 
95 and /n TO , s and /. i Wt i s , their chemical potentials. Equation (2) can be reorganized as, 

G 2 = X w ^l w 2 b“ -^yMy, 2 “b IT's (M w,s Mtu,2) “b H Is (M w,ls M^, 2 ) (3) 

Using Eqs. (1) and (3) the work of germ formation, AG = G 2 — G\, can be written as 

AG — AG s l n ~\~Tl s (^fX W:S k’w, 2 ) “b Tll s ([A w ,ls Mu;, 2 ) (4) 

where AG s i n is the change in the Gibbs free energy of the bulk solution caused by the appearance 
of the germ, i.e., 

AG s l n ~ Au;(Mu;,2 Mu;, l) “b A^y(My,2 My, l) (3) 

Eq. (4) indicates that the work of germ formation originates from (i) changes in the composition of 
100 the liquid phase, (ii) the formation of the interface and (iii) the formation of the bulk of the solid. 
Using the equilibrium between ice and the liquid solution as reference state, the latter can be written 
in the form (Kashchiev, 2000), 

Vw, s ~ Vw,2 = -k B T\n[ ( 6 ) 

\ ^w : eq J 

where k B is the Boltzmann constant and a Wt< , q is the equilibrium water activity between bulk liquid 
and ice. 

105 AG s in in Eq. (5) arises because the solute must be “unmixed” (Black, 2007) from the liquid to 
form a solute-free germ, causing a change in the molar composition of the liquid phase. Such unmix- 
ing carries an entropic cost to the system (Bourne and Davey, 1976). Thus, AG s i n is proportional to 
the mixing entropy of the atoms in the germ, 

AG s in = nk B T (X w lnUn, Aylndy) (7) 

Where X w and X y are the molar fractions of water and solute in the bulk liquid, respectively. For 
1 1 0 dilute solutions, which is almost always the case for ice nucleation, X w » X y and, 

AG s \ n ^ -nk B T\na w (8) 

Where n = n s + ni s is the total number of molecules in the germ. 
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2.1 Energy of formation of the interface 


To further develop Eq. (4) it is necessary to introduce a model of the solid-liquid interface. Theo- 
retical models show that the solid-liquid interface is characterized by the organization of randomly 
moving liquid molecules into positions determined by the solid matrix (Spaepen, 1975; Karim and 
Haymet, 1988; Haymet and Oxtoby, 1981). Associated with this increased order is a decrease in 
the partial molar entropy of the liquid molecules. Since the solid determines the positions of the 
molecules at the interface, the partial molar entropy at the interface must approximate the bulk en- 
tropy of the solid. However the interface molecules are liquid-like, and their enthalpy remains that 
of the bulk liquid (Black, 2007). This implies that the system must pay the maximum entropic cost 
during the formation of the germ (Spaepen, 1975; Black, 2007). The entropic nature of the thermo- 
dynamic barrier for nucleation has been confirmed by molecular dynamics simulations (Reinhardt 
and Doye, 2013). This conceptual model is used below to develop an expression for the energy of 
formation of the interface. 

The change in the partial molar free energy of water associated with the formation of the interface 
is given by 


w,ls — h w J s Ts w J s 2 (9) 

Where s w j s is the entropy of the interface molecules. Assuming that the entropy of the molecules at 
the interface approximates the entropy of the bulk solid, i.e., s w j s ~ s w>s , Eq. (9) can be written as, 

w : ls l^w : 2 hw,ls Ts w ,s fJ"w,2 (10) 

Taking into account that fi W:S = h WtS —Ts w , s , and using Eq. (6) into Eq. (10) we obtain 

ft w,is ^w,2 — /cgTln ( j + \h Wi i s (11) 

where A h w ^ s = h Wt i s — h WiS is the excess enthalpy of the molecules at the interface. 

If no solute is present the enthalpy of the molecules at the interface approximates the enthalpy of 
water in the bulk, i.e, A h w j s « A hf where A hf is the latent heat of fusion of water. However the 
adsorption of solute at the interface affects Ah w f s . Using the Gibbs model of adsorption the effect 
of the solute on A h w j s can be accounted for in the form (Hiemenz and Rajagopalan, 1997; Gibbs 
et al„ 1928), 


A h w f s = A h f - V w k^T\na w - T y kBTlna y (12) 

where P,,. and T y are the surface excess of water and solute, respectively, and represent the ratio 
of the number of molecules in the interface to the number of molecules at the dividing surface. 
According to the Gibbs model, T,,, and T y depend on the position of the dividing surface, typically 
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chosen so that the surface excess of solvent is zero (Kashchiev, 2000). However by choosing the 
140 dividing surface as equimolecular with respect to the solute (Section 2) the resulting expressions 
become independent of the nature of the solute. Thus making T y = 0, Eq. (12) becomes, 

Ah w> u = Ahf-T w k B T\n.a w (13) 

In the solid matrix the number of molecules at the surface is given by sn 2 / 3 where s is a geometric 
constant depending on the crystal lattice (1.12 for hep crystals and 1.09 for bcc crystals (Jian et al., 
2002)), and n is the total number of atoms in the germ. However the interface is generally made of 
145 several layers beyond the outer layer of the solid (Henson and Robinson, 2004; Chen and Crutzen, 
1994; Spaepen, 1975). Spaepen (1975) showed that for random coverage of a solid matrix there are 
about 1.46 molecules at the interface for each molecule in the outer layer of the solid matrix. With 
this, r„, = 1.46s and ni s = 1.46sn 2 ' 3 . Equation (13) then becomes. 


Ah w j s = Ahf — lA6sk B Tlna w (14) 

Introducing Eq. (14) into Eq. (1 1) we obtain, 

Hw,ls ~ Hw,2 = -k B T\n( -°A— J +A/i / -1.46sfc B 2 n lna u; (15) 

\ ®’ w,eq J 

150 Equation (15) expresses the energy cost associated with the formation of the interface accounting 
for solute effects. Since it results from the consideration of the maximum entropy increase across 
interface, this model will be referred to as the Neg-Entropic Nucleation Theory (NENT). 

2.2 Nucleation work 


Introducing Eqs. (6), (8) and (15) into Eq. (4), and rearranging we obtain 

A G = -nk B T\n ( + 1.46sn 2/3 (A h f - lA6sk B T\na w ) 

\ Q"iv,eq J 


(16) 


1 55 where n = n s + nu was used. 

The germ size at nucleation, n* , and the nucleation work, AG nuc , are obtained by applying the 
condition of mechanical equilibrium to Eq. (16), i.e.. 


dAG nuc 

dra* 


-feBTln^ — — ^ +-1.46s(n*) x ^ 3 {Ahf— H.46sfcBTlna tl; ) = 0 (17) 

V ®w,eq J o 


Solving Eq. (17) for n* and rearranging gives. 
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1.46s(A/i / — l.TbsfcByinauj) 


(18) 


The nucleation work is obtained by replacing Eq. (18) into Eq. (16). After rearranging we obtain, 
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AG nU c — 


4 [1.46s(Ahf — lAGskBThian,)] 3 


27 


1 60 The nucleation rate, Jhom, is given by. 




(19) 


Jhom = exp ^ - A ^ C j (20) 

where Jo is a T-dependent preexponential factor. As in CNT, it is assumed that Jo results from 
the kinetics of aggregation of single water molecules to the ice germ from an equilibrium cluster 
population (Kashchiev, 2000), therefore. 


Jo = 


N c k B T p w ZQ g 


-exp - 


A G a 


( 21 ) 


h Pi v w V k ^T 

where N c is the number of atoms in contact with the ice germ, p w and pi are the bulk liquid water 
1 65 and ice density, respectively, is the germ surface area, and AG act is the activation energy of the 
water molecules in the bulk of the liquid. Z is the Zeldovich factor, given by 


Z = 


AGn 


3nk B Tn* 


1/2 


( 22 ) 


2.3 Classical Nucleation Theory 


CNT is commonly used to described homogeneous ice nucleation (Khvorostyanov and Curry, 2004) 
and is therefore important to compare the NENT model against CNT predictions. According to 
1 70 CNT, the nucleation rate is given by (Pruppacher and Klett, 1997) , 


^CNT 


2 N c k B T p w 

h pi 


( &iw \ 1/2 

AGact + AGcnt 

exp 

k B T 


The work of nucleation, AGcnt, is calculated as 


(23) 


AGcnt = 


167 raf, 


(24) 


3(fc B Tln5 l ) 2 

where Si = a w (p s ,w/Ps,i)> is the saturation ratio with respect to the ice phase. The critical germ size 
is given by, 


n CNT — 


(25) 


3(k B TlnS t ) 3 

The usage Eq. (23) requires a parameterization of (T; w , for which there is a large uncertainty. The- 
1 75 oretical approaches have been developed to estimate cr; w however they are mostly applicable to low 
undercooling (e.g., Digilov, 2004; Spaepen, 1994) and (j; w is in general found by fitting experimen- 
tal measurements of Jhom (e.g., Murray et al., 2010a; Khvorostyanov and Curry, 2004; Marcolli 
et al., 2007). Two approaches are employed to parameterize <7 iw . Following Murray et al. (2010a), 
the following correlation was used to describe <t; w , 
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(T) = 0.0229 f^—) (Jm“ 2 ) 

v ’ V 236.8 J v ’ 


( 26 ) 


With T in K. The parameters of the correlation in Eq. (26) were slightly modified to match the 
freezing point of pure water measured by Koop et al. (2000). The model presented in Section 2 
however indicates that besides T, er iw must also depend on a w since the presence of the solute in the 
interface layer modifies the interfacial energy. To account for this, a correlation for er; w was obtained 
by fitting Eq. (23) to the data of Koop et al. (2000) in the form. 


a iv ,(T,a w ) = 0.00211 -0.0513a™ + 3.04 x 10“ 4 T (J m -2 ) (27) 

With 180 <T < 273 in K. The linear dependency of cri w on T is consistent with theoretical stud- 
ies (Spaepen, 1994; Pruppacher and Klett, 1997). In agreement with experimental measurements 
(Pruppacher and Klett, 1997; Digilov, 2004), Eq. (27) predicts o-; w = 33.9 J m -2 for T = 273 K and 


3 Discussion 
3.1 Nucleation Rate 

Figure 2 shows the nucleation rate calculated from the NENT and CNT models, i.e., Eqs. (20) 
and (23). For NENT, the surface area parameter, s, in Eq. (19) was set to 1.105, that is, the ice 
germ structure is assumed to lie somewhere between a bcc (s = 1.12) and a hep (s = 1.109) crystal, 
justified on experimental studies showing that ice forms as a stacked disordered structure (Malkin 
et al., 2012). The values used for the parameters of Eqs. (19) and (20) are listed in Table 1. Jh om 
from the empirical correlations of Koop et al. (2000) and Murray et al. (2010a) is also shown in Fig. 
2. Murray et al. (2010a) compared experimentally determined nucleation rates from several sources 
and found about a factor of 10 variation in Jhom of pure water. The correlation of Murray et al. 
(2010a) was therefore included in Fig. 2 for reference, although it is only applicable around 236 K. 
CNT results are shown using Eqs. (26) and (27) to calculate <r iw . 

The median freezing temperature, T me( j, was calculated by solving, 


Jhom(Tjned)AiUd = 1 (28) 

where At is the experimental time scale va the droplet volume. Equation (28) is a result of the 
nucleation dispersion theory (Barahona, 2012) and indicates that half of the droplets in a population 
will freeze when there is on average one ice germ per droplet produced during At. T me d is calculated 
by numerical iteration and assuming At = 10 s and a droplet diameter of 10 /im, selected to match 
to the conditions used by Koop et al. (2000). 
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There is overlap between all the curves of Fig. 2 for T around 236 K, that is, near the homogeneous 
freezing T of pure water ( a w = 1). This is in agreement with the study of Murray et al. (2010a) 
showing that most models predict similar Jh om for pure water. However for T < 236 K and for 
a w < 1, Fig. 2 shows significant differences in predicted nucleation rates. At constant a w , Jhom 
for NENT increases more rapidly as T decreases than for CNT, indicating a more rapid decrease in 
AG nuc in NENT. Jh om for NENT also increases more quickly than for K00 as T decreases, although 
the discrepancy remains within the typical scatter of experimentally determined Jhom (Murray et al., 
2010b). When Eq. (27) is used to parameterize cri w , Jhom from CNT shows good agreement with 
K00; this is however by design as K00 was used to generate Eq. (27). When Eq. (26) is used to 
parameterize a lvl , Jhom from CNT is much lower than predicted by either NENT or K00, and only 
at a w = 1 Jhom from CNT agrees with experimental observations. Such high sensitivity to er iw is 
one of the main drawbacks of CNT. 

At constant a w , NENT and CNT show an initial increase in Jhom as T decreases however this ten- 
dency reverses at low T. This behavior is caused by an increase in AG act as T decreases, reducing 
Jhom- As T decreases the role of activation of water molecules becomes increasingly more signifi- 
cant, limiting Jhom- In contrast, the K00 parameterization shows a monotonically increase in Jhom 
as T decreases. A plausible explanation for this behavior may be found in the experimental method 
used to generate the K00 parameterization. Optical methods are accurate at low Jhom where the 
freezing of individual droplets can be easily discriminated. However as T decreases, high Jhom re- 
sults in rapid freezing of the entire droplet population. Thus the maximum and subsequent reduction 
in Jhom as T decreases may be difficult to infer from observation of the freezing of single droplets. 
Thus it is likely that the K00 parameterization overestimates the highest values of Jhom- Numerical 
test (not shown) suggest this may lead to overestimation of droplet freezing fractions particularly for 
small droplets, although it may have a limited impact in ambient clouds where freezing fractions are 
typically small (Barahona and Nenes, 2008). 

3.2 Critical germ size 

Figure 3 shows the critical germ size in terms of number of water molecules, calculated using NENT 
and CNT. According to the nucleation theorem (Kashchiev, 2000), n* can also be determined di- 
rectly from experimental results following, 

* dAG nuc 

Ti 

dA/j, w d/\[i w 

where $ is the energy of formation of the interface, and A fj, w 
can be simplified as (Kashchiev, 2000), 


(29) 

= — kBT\n(^ a aw ^ . Equation (29) 


n 


* 


din Jhom 
dlno^ 


(30) 
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where it is assumed that $, or equivalently cr iw , does not depend on a w . Using the K00 parameter- 
ization into Eq. (30) results in n* between 400 and 600 molecules for T between 180 K and 236 
K (Fig. 3). On the other hand, NENT (Eq. 18) predicts n* around 300 for the same T interval 
and CNT around 100 (Eq. 25). A similar discrepancy was found by Ford (2001) who ascribed it to 
deficiencies in CNT. However NENT offers further insight into the origin of the differences in n*. 

Figure 4 suggests that around the median freezing a w (defined in a similar way as in Eq. 28), 
d eUna° m f° r K00 is similar as for CNT (with <r iw as defined in Eq. 27) and lower than for NENT. The 
fact that n* is higher for K00 than for NENT (Fig. 3), even though is higher in the latter, is 

at odds with the predictions of Eq. (30). This picture can only be reconciled if ^ 0, that is, the 
interfacial energy must be a functional of a w , which is also suggested by Eq. (19). This dependency 
is explicit in the NENT model (Section 2.1). In CNT, it is introduced by making erj w a function of 

Qj w . 

To explain the dependency of the interfacial energy on a w one must consider the Gibbs model 
of the interface (Section 2.1). By introducing the arbitrary dividing surface, an excess number 
of molecules is created around the interface between the liquid and the solid (Hiemenz and Ra- 
jagopalan, 1997). This is typically dealt with by selecting the so-called equimolecular dividing 
surface (EDS), in which the interface has energy but its net molecular excess is zero (Kashchiev, 
2000; Schay, 1976). However the EDS cannot be defined simultaneously for both, the solute and the 
solvent. In fact, using the EDS with respect to the solvent, as implicitly done in CNT, results in a 
molecular excess of solute at the interface. This explains the discrepancy in n* between Eq. (30) and 
CNT. In Section 2 it is shown that it is advantageous to define the EDS with respect to the solute, and 
account explicitly for the excess of water molecules at the interface. Thus the consistency between 
the choice of the dividing surface and the molecular excess at the interface is explicit in NENT. This 
picture also implies that Eq. (29) (which does not depend on the choice of the dividing surface) 
instead of Eq. (30), must be used in the analysis of ice nucleation data. 

3.3 Median Freezing Temperature 

Finally we investigate whether the model presented in Section 2 is able to explain the water activity 
criteria of Koop et al. (2000). Figure 5 shows the median freezing T defined by Eq. (28), calculated 
using K00, CNT and NENT. For CNT using Eq. (27) results in agreement with K00, which however 
is by design as T me( j predicted by K00 was used to specify (j; w (Eq. 27). Using er iw from Eq. (26) 
which is based on a fit to observed T me d at a w = 1 results in overprediction of T me d for a w < 1. 
T me a from K00 and NENT overlap down to 190 K (Fig. 5). Since no data from K00 are used in 
NENT, the model developed here constitutes the first theoretical explanation of the results of Koop 
et al. (2000). 

The shift in water activity, A a w =a w — a Wieq , can be obtained by solving, 
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fc B Tln(J 0 Ato d ) 


(31) 


4 {1.46s[A/iy — lAGsk-QT\w(a w ^ q -\- Aa^)]}' 3 
27 {fc B Tln[ ^+f ;^ 2 ]} 2 

Equation (31) was obtained by replacing Eq. (20)intoEq. (28). Since the roots ofEq. (31) determine 
T me d, it is termed the characteristic freezing function . 

Inspection of Eq. (31) shows that it is a function of T only, since the a w dependency is removed 
by application of Eq. (28). Thus, the roots of Eq. (31) are determined by the value of Aa^. Figure 
(6) shows that Eq. (31) only has real solutions in the interval 185 < T < 238 K over a very narrow 
set of values of A a w , i.e., 0.298 < A a w < 0.304. This is the origin of the water activity criteria since 
such variation in A a w is well within experimental uncertainty. An interesting feature of Eq. (31) is 
that it produces similar T — a w curves for different A a w values. This means that the multiple roots 
of Eq. (31) are located at similar T med for different values of Aa w , and always fall on the same 
curve. 

Figure 6 shows that Eq. (31) constitutes a theoretical derivation of the water activity criteria. Aa„, 
can be easily obtained by numerical iteration of Eq. (31). However at a w = 1, Eq. (31) is greatly 
simplified and A a w can be found by direct analytical solution, in the form. 


Aa w = 1 — exp 


2 / Ahf \ 3/2 

3 x /31n(J 0 Afn d ) V T* ) 


= 0.2981 


(32) 


where T* = 236.03 is the median freezing temperature at a w = 1. The value of Aa w in Eq. (32) was 
obtained using the parameters of Table 1 calculated at T* . A a w is very close to the experimental 
value of 0.302 found by Koop et al. (2000). For T > 190 K, A a w calculated from Eq. (31) is fairly 
constant (being 0.300 at T = 190 K). For T < 190 there is a slight increase in A a w reaching about 
0.31 at T = 180 K. This increase is due to the increase in AG act at low T. Jhom at very low T is 
still uncertain, since factors like the formation of glasses (Murray et al., 2010b) and the formation 
of highly concentrated brines within droplets may play a role (Bogdan and Molina, 2010). The 
accuracy of K00 at very low T has also been questioned (Swanson, 2009). 

The preexpoential factor in Eq. (21) is almost constant between 180 K and 236 K. That is, the 
flux of molecules to the germ is controlled by A G act . In fact, by replacing Eq. (19) into Eq. (20) 
and then into Eq. (28), we obtain after rearranging. 


AG nuc + AG ac t 
T 


constant 


(33) 


Thus an increase AG act is balanced by a decrease in AG nuc , i.e., the increase in the driving force 
for nucleation at low T makes up for the decrease in the mobility of water molecules. Since AG nuc 
can be defined over a pure thermodynamic basis (Section 2), Eq. (33) suggest that AG act and the 
kinetic term in Eq. (20) may also admit a thermodynamic description. 
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4 Conclusions 


The model presented in Section 2 constitutes a new theory of nucleation that does not use the in- 
terfacial tension as defining parameter. It is therefore free from biases induced by uncertainties in 
the parameterization of <r iw . Instead, an expression for the interfacial energy was developed directly 
using thermodynamic principles. The new theory is based on a conceptual model in which the in- 
terface is considered to be made of “water molecules trapped by the solid matrix'’. It also accounts 
for the finite droplet size leading to changes in the composition of the liquid phase upon nucleation. 
Since it places emphasis on the increase in order and the reduction in entropy across the interface, 
the new model has been termed the Neg-Entropic Theory of Nucleation, NENT. 

Comparison against experimental results showed that the new theory is able to reproduce mea- 
sured nucleation rates and is capable of explaining the observed constant shift in water activity 
between melting and nucleation (Koop et al., 2000). The constant water activity shift originates be- 
cause the median freezing temperature only exist for a very narrow range of A a w values (Eq. 31). 
More fundamentally, NENT shows that the effect of water activity on nucleation is a manifestation 
of the entropic bander to the formation of the germ. An analytical expression for Aa„, was derived 
and was shown to agree well with the experimental value measured by Koop et al. (2000). This 
constitutes the first theoretical derivation of the water activity criteria found by Koop et al. (2000). 

The new theory suggests that to reconcile experimental results and theoretical models the interfa- 
cial energy must depend on a w . This is implicit in the development of NENT, however is missing 
in CNT. The dependency of Jhom on a w originates from the excess concentration of either solute of 
solvent when the dividing surface is defined. Such excess is present even if the EDS is defined with 
respect to the solvent. It was shown that it is advantageous to define the EDS with respect to the 
solute because the resulting expressions are independent of the nature of the solute, and therefore 
consistent with experimental observations. Although such considerations are neglected in CNT, it 
was shown that CNT can be corrected by allowing cr iw to depend on a w , for which a new empirical 
correlation was developed. 

Analysis of the new theory suggested that the temperature dependency of both the kinetic and 
thermodynamic terms plays a significant role in defining Jhom and T me d- It was shown that around 
Jmed the increase in AG ac t as T decreases is compensated by a decrease in AG nuc . Thus an 
increased driving force for nucleation compensates the slower molecular diffusion at low T. Such 
coupling between kinetic and thermodynamics during nucleation suggests that a thermodynamic 
description of the preexponential factor (Eq. 20) may be possible. 

Some disparity was found between the K00 parameterization and NENT for T < 190 K. This may 
be ascribed to non-equilbrium effects and the possible formation of glasses at low T. Also the solute 
term in Eq. (7) may become important as a w decreases. Further research is required to elucidate 
the mechanism of freezing at such low T. The model presented here emphasizes the entropic nature 
of homogeneous nucleation. Molecular simulations may shed further light on the role of entropy 
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changes across the interface on ice nucleation. Measurements of the interface thickness would also 
help elucidate the role of the ice crystal geometry (represented by the constant s in Eq. 15) in 
340 determining Jhom- The new theory is also suitable to describe heterogeneous freezing, particularly 
in the immersion mode. This will be addressed in future studies. 

The theory presented here reconciles theoretical an experimental results. It will also help reduce 
the uncertainty in Jhom associated with the parameterization of cr; w in theoretical models. It offers 
for the first time a thermodynamically consistent explanation of the effect of water activity on ice 
345 nucleation. Its relative simplicity makes it suitable to describe ice nucleation in the atmosphere, 
and may lead to a better understanding of the formation of ice clouds and their parameterization in 
atmospheric models. 
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Table 1 . List of symbols. 


CLw ? CLy 

Activity of water and solute, respectively 

f^w,eq 

G 

Equilibrium a w between bulk liquid and ice (Koop et al., 2000) 
Gibbs free energy 

h 

Planck’s constant 


h w , s , h W 'i s Partial molar enthalpy of water in ice and at the interface, respectively 
Jo Preexponential factor 


J hom 

Nucleation rate 

&B 

Boltzmann constant 

n 

Total number of molecules in the solid germ 

n* 

Critical germ size 

Tils 

Number of molecules in the bulk of the solid and in the interface, respectively 

N c 

Number of atoms in contact with the ice germ, 5.85 x 10 18 m~ 2 (Pruppacher and Klett, 1997) 

N W ,Ny 

Total number of water and solute molecules, respectively 

Ps,w, Ps,i 

Liquid water and ice saturation (Murphy and Koop, 2005) 

s 

Geometric constant relating n and ni s 

Si 

Saturation ratio with respect to ice 

Sw,s 5 Svu,ls 

T 

Partial molar entropy of water in bulk ice and at the interface, respectively 
Temperature 

Tined 

Median freezing T 

Vd 

Droplet volume 

Vw 

Molecular volume of water in ice (Zobrist et al., 2007) 

Xyj , Xy 

z 

Molar fraction of water and solute, respectively 
Zeldovich factor 

AG ac t 

Activation energy of liquid water (Zobrist et al., 2007) 

AGnuc 

Nucleation work 

AGsln 

Change in free energy of the bulk solution during nucleation 

A h W 'i a 
A hf 
At 

Excess enthalpy of the interface 
Heat of fusion of water a 
Experimental time scale 

Adw 

Giui &w,eq 

<s> 

Energy of formation of the interface 

r„, r y 

Molecular surface excess of water and solute, respectively 

f-lw 5 l^y 

Chemical potential of water and solute, respectively 

f^w,ls 

Chemical potential of water at the interface 

l-lw,S 

Chemical potential of bulk ice 

pwi Pi 

Bulk density of liquid water and ice, respectively (Pruppacher and Klett, 1997) 

(Jiw 

Ice-liquid interfacial energy 

Clg 

Ice germ surface area 


"Fmrn the data of Johari et al. (1994) the following fit was obtained: Ahf — 7.50856 X 10 7 T 5 — 
8.40025 X 10 -4 T 4 + 0.367171X 3 — 78.1467T 2 + 8117. 02T — 3.29032 X 10 s (J mol -1 ) for T be- 


tween 180 K and 273 K. 
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Fig. 1 . Scheme of the formation of an ice germ from a liquid phase. 
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Fig. 2. Homogeneous nucleation rate. K00 and M10 refer to results obtained using the correlations of Koop 
et al. (2000) and Murray et al. (2010a), respectively. For CNT <Ti W was calculated using either Eq. (27) (line 
CNTa) or Eq. (26) (line CNTb). Lines are grouped by increasing water activity: a w = 0.8, 0.9 and 1.0, from 
left to right, respectively. 
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Fig. 3. Critical germ size, n*. 
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Fig. 4. Homogeneous nucleation rate. The freezing line was calculated using Eq. (28). Lines are grouped by 
temperature: T = 195,220 and 236 K, from left to right, respectively. 
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Fig. 5. Median freezing temperature for homogeneous nucleation. 
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Fig. 6. Characteristic freezing function, Eq. (31). 
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